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1.0  INTRODUCTION 


\ 


—  This  presentation  is  a  summary  of  digital  ionogram 
processing  techniques. developed  by  ULCAR  in  the  recent  past. 
Real-time  and  on-line  methods  were  emphasized  and  this  implies 
continuing  improvement  and  increasing  complexity  of  the  proc¬ 
essing  algorithms  along  with  the  technological  advancements. 

Of  course,  a  mini  or  mainframe  computer  with  sufficient  memory 
could  perform  the  processing  described  here.  But  for  on-line 
operation  at  an  ionosonde  field  site  cost  and  environmental 
conditions  generally  exclude  this  approach,  and  special  pur¬ 
pose  hardware/microprocessor  techniques  are  the  only  adequate 
method.  The  difference  between  mini-  and  microcomputers  is 
diminishing,  but  the  modular  microcomputer  systems  have  the 
advantage  of  being  easily  integrated  into  a  digital  ionosonde 
at  low  costs .  The  Digisjonde  128PS.  (Bibl  and  Reinisch,  1978a; 
1978b)  initially  outputs  complete  ionograms  (in  the  ionogram 
mode)  with  amplitude,  phase,  doppler,  incidence  angle  and/or 
polarization  for  each  frequency-range  bin.  Post-processing  of 
these  ionograms  can  then  extract  the  ionospheric  parameters 
required  for  specific  investigations. 


2.0  AUTOMATIC  SCALING  OF  IONOGRAMS 


The  initial  ionogram  contains  a  large  amount  of  in¬ 
formation  of  which  not  all  is  necessarily  required  for  syn¬ 
optic  studies  or  investigations  interested  only  in  specific 
parameters.  In  such  cases  it  may  suffice  to  "scale”  the  iono¬ 
gram  by  automatically  determining  the  heights  and  amplitudes 
of  the  main  echo  traces.  For  undisturbed  ionograms  the  task 
is  relatively  simple  as  long  as  the  distinction  between  or¬ 
dinary  and  extraordinary  echoes  does  not  become  too  difficult. 
Figure  1  shows  the  result  of  a  simple  algorithm  that  was  ap¬ 
plied  to  a  mid-latitude  ionogram.  This  algorithm,  implemented 
in  the  Geomonitor  (Reinisch  and  Smith,  1976)  that  is  connected 
to  the  Digisonde  at  AFGL’s  Goose  Bay  Ionospheric  Observatory, 
Labrador  (Buchau  et  al,  1978)  and  in  the  Automatic  Ionogram 
Collator  at  the  U.S.  Army  Digisonde  in  Fort  Monmouth,  New 
Jersey,  extracts  amplitude,  height  and  spread  of  the  echoes 
from  the  ionogram.  These  collated  data  can  be  presented  in 
the  frame  of  an  ionogram  to  verify  the  accuracy  of  the  algo¬ 
rithm,  as  shown  on  top  of  Figure  1,  or  they  can  be  printed  in 
the  form  of  ionospheric  characteristics.  Figure  2  shows  E 
and  T  echo  amplitudes  as  function  of  frequency  and  time  (panels 
2,  3,  4  and  5),  and  the  range  for  backscatter  observations 
(panel  1)  as  function  of  time  for  a  period  of  s.ix  days  in 
Goose  Bay.  Each  vertical  line  contains  the  data  of  one  iono¬ 
gram.  The  envelopes  of  the  amplitude-versus-frequency  char¬ 
acteristic  give  essentially  the  diurnal  variations  of  the 
critical  frequencies  for  E  and  F-region.  Disturbed  iono¬ 
spheric  conditions,  like  on  days  291/292  1977  are  easily 
recognizable. 

Most  ionograms  are,  of  course,  not  as  simple  as  the 
one  shown  for  a  mid-latitude  station  in  Figure  1.  Figure  3 
shows  two  somewhat  more  complicated  ionograms  from  AFGL's 
Goose  Bay  Ionospheric  Observatory  and  the  "echo"  points  de¬ 
tected  by  the  on-line  Geomonitor.  Of  special  interest  is 


the  upper  ionogram  from  10  January  1979  which  contains  a 
strong  Z-trace.  Now,  even  though  the  new  Digisonde  128PS  tags 
each  frequency-range  bin  with  an  0  or  X-flag  the  Z-trace  is 
not  distinguished  from  the  0-trace  since  both  have  the  same 
polarization.  Since  the  Geomonitor  detects  up  to  six  frequen¬ 
cies  for  each  sounding  frequency  it  will  pick  up  both  the  0 
and  the  Z-trace  and  the  Automatic  Parameter  Evaluation  (A.P.E.) 
program  (Section  3.0)  will  determine  a  false  foF2.  In  general 
further  processing  is  required  to  extract  a  unique  trace  that 
can  serve  as  input  to  the  true  height  algorithm.  It  is  clear 
that  a  refined  ionogram  (or  the  corresponding  profile)  has 
lost  information  contained  in  the  original  ionogram,  like 
oblique  echoes  and  spread  conditions.  It  is  for  that  reason 
that  for  many  investigations  reliance  on  a  few  ionospheric 
parameters  is  dangerous,  and  the  raw  ionogram  data  must  be 
preserved.  On  the  other  hand  for  routine  monitoring  of  iono¬ 
spheric  conditions  limitation  to  characteristic  parameters 
becomes  a  necessity. 
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3.0  AUTOMATIC  PARAMETER  EVALUATION 


Once  the  echo  traces  are  determined  it  is  possible 
to  determine  certain  ionospheric  parameters  which  in  conven¬ 
tional  systems  are  manually  scaled.  The  Automatic  Parameter 
Evaluation  program  (Smith  et  al,  1979)  was  applied  to  some 
1000  Geomonitor  processed  Goose  Bay  ionograms  to  extract  the 
parameters  foF2,  fmin,  hmin,  MUFC3000)  and  ftEs.  For  the 
period  5  to  13  January  1979  the  automatically  scaled  values 
were  compared  with  a  data  technician’s  visual  interpretations. 
Figures  4  and  5  show  the  diurnal  variations  of  the  errors  for 
foF2  and  the  MUFC3000).  Considering  that  spread  F  occurrence 
and  mid-latitude  trough  phenomena  effected  the  majority  of 
the  ionograms  during  the  5  to  13  January  period,  the  success 
rate  shown  in  Figures  4  and  5  is  fairly  good  during  daytime 
and  around  local  midnight.  Performance  during  twilight  times 
is  unacceptably  poor.  The  weak  point  in  the  automatic  proc¬ 
essing  scheme  is  the  presently  used  Geomonitor  algorithm  for 
the  on-line  echo  trace  detection  under  disturbed  conditions. 

An  improved  algorithm  applicable  to  both  bottomside  and  top¬ 
side  ionograms  (Reinisch  et  al,  1980)  is  presently  being  de¬ 
veloped  using  the  ionograms  as  an  entity  rather  than  making  a 
frequency-by-frequency  decision  on  what  constitutes  an  echo. 
The  recent  advancements  in  memory  technology  allow  on-line 
microprocessor  implementation  of  this  advanced  algorithm. 

Automatic  determination  of  the  predicted  maximum 
usable  frequency  (MUF) ,  for  a  given  distance,  from  the  verti¬ 
cal  ionograms  is  not  difficult  for  digital  ionograms.  Methods 
to  determine  MUF  as  function  of  time  had  been  developed  for 
analog  sounders  (Nakata  et  al,  1953;  Bibl,  1956).  Our  digital 
algorithm  (Smith  et  al,  1978)  simply  multiplies  each  frequency 
of  the  ionogram  trace  by  the  MUF  factor  calculated  for  a  dis¬ 
tance  of  3000  km  according  to  URSI  specifications  (Piggott 
and  Rawer,  1972,  p.  21)  to  generate  the  equivalent  oblique 


ionogram  (Figure  6)  from  which  the  MUF  can  be  scaled  directly. 
More  recently  we  applied  this  method  to  the  raw  ionograms 
rather  than  to  the  Geomonitor  identified  traces.  The  average 
error  diminished  because  errors  in  the  trace  identifications 
no  longer  effects  the  MUF  determination.  Cormier  and  Dieter 
(1974)  simulated  the  manual  overlay  technique  with  good  suc¬ 
cess  in  a  computer  program  to  find  the  MUF.  This  method  is, 
however,  more  susceptible  to  noise  than  the  approach  discussed 
above . 


4.0  REAL  TIME  ELECTRON  DENSITY  PROFILES 


For  conversion  of  the  h'(f)  ionogram  trace  into  ver¬ 
tical  electron  density  profiles  the  Digisonde  contains  a  micro¬ 
computer  program  that  performs  the  conversion  in  almost  real 
time.  The  algorithm  is  based  on  the  usual  lamination  techni¬ 
que  for  a  monotonous  profile  assuming  that  z(fn>,  i.e.  height 
as  function  of  plasma  frequency,  is  continuous  and  steady  at 
slab  transitions  and  that  the  second  derivity  of  z  with  re¬ 
spect  to  fN  is  constant  within  each  slab  (Paul,  1977;  Bibl 
and  Reinisch,  1976b,  p.  98ff).  This  results  in  a  recursive 
scheme  that  determines  the  height  of  plasma  frequency  f  as: 


z  =  z  ,  +  D  ,+i-DD  , 
n  n-1  n-1  2  n-1 
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The  group  index  of  refraction  y*(f,f  )  is  of  course  also  a 

function  of  the  gyrofrequency  f^.  For  the  on-line  profile 

calculations  the  coefficients  F  .  and  S  .  are  tabulated  for  a 

ni  ni 

given  geomagnetic  location  and  stored  in  PROM  to  achieve  a 
fast  profile  conversion  in  the  Microcomputer. 


At  the  present  time  the  ionogram-to-prof ile  conver¬ 
sion  is  not  fully  automated.  The  f  ,  h*  data  pairs  must  be 

n  n 
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entered  via  keyboard.  An  example  of  the  program  input  and 
output  data  is  given  in  Table  1.  These  on-line  electron  den¬ 
sity  profiles  are  of  great  importance  for  special  experiments 
like  rocket  campaigns.  The  example  shown  was  made  during  a 
rocket/satellite/incoherent  radar  experiment  in  Kwajalein  on 
31  July  1978.  The  last  column  in  the  table  shows  the  bottom- 
side  scale  length  of  the  F-layer,  giving  the  height  interval 
within  which  the  plasma  frequency  increases  by  1  MHz.  This 
scale  height  together  with  hminF  seems  to  control  the  spread  F 
occurrence  (Ossakow  et  al,  1978).  No  effort  was  made  to  re¬ 
fine  the  program  with  regard  to  the  starting  and  valley  prob¬ 
lems.  Instead  we  optimized  for  speed  and  ease  of  operation  at 
the  field  site. 
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INPUT  DATA,  KW  78  0212  2300  LMT  TRUE  HEIGHT  PROFILE 

[R  FREQ.  [100  kHz]  V.  HEIGHT  [KM]  Z  [KM]  SCALE  [KM/MHz] 
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TABLE  1 


5.0  DOPPLER  IONOGRAMS 


The  usefulness  of  multiparameter  ionograms  (Bibl 
and  Reinisch,  1978a)  is  demonstrated  by  the  ionograms  re¬ 
corded  by  AFGL's  Airborne  Ionosonde  in  Thule  on  9  December 
1979  (Figure  7).  The  amplitude  ionograms  in  the  lower  right 
hand  corner  shows  typical  aurora  region  spread  F  conditions. 
The  amplitude  in  each  frequency-range  bin  is  accompanied  by  a 
status  word  displayed  in  the  upper  right  corner,  forming  a 
status  ionogram.  In  this  instance  the  status  words  represent 
signal  doppler  frequencies  varying  from  -11  to  +11  Hz.  A 
simple  microcomputer  algorithm  separated  the  original  ionogram 
into  two:  one  containing  only  echoes  with  negative  doppler 
frequencies  (upper  left),  and  another  with  only  positive  dop¬ 
pler  frequencies  (lower  left).  In  these  doppler  ionograms 
the  amplitudes  are  replaced  by  doppler  frequencies.  The  over¬ 
head  ionosphere  shows  no  significant  motion,  slightly  positive 
and  negative  dopplers  are  observed.  Predominantly  negative 
doppler  frequencies  are  observed  for  ranges  between  400  and 
600  km.  The  radial  velocity  of  the  reflecting  ionization  away 
from  the  sounder  is  110  ms”1.  Between  8  and  9  MHz  echoes 
with  a  strong  positive  doppler  shift  are  seen  at  ranges  be¬ 
tween  600  and  700  km.  The  doppler  frequencies  indicate  a 
radial  velocity  of  the  reflector  of  approximately  75  ms”1. 

If  the  status  contains  information  on  the  incidence 
angle  of  the  ionospheric  echoes  directional  ionograms  can  re¬ 
veal  the  structure  of  the  local  ionosphere.  Figure  8  shows 
two  directional  ionograms  from  the  1978  Kwajalein  equatorial 
spread  F  campaign.  The  receiving  antenna  array  was  scanning 
the  vertical,  northwest  and  southwest  direction  in  the  first 
ionogram,  and  vertical,  northeast  and  southeast  in  the  second 
one.  The  majority  of  the  echoes  come  from  overhead,  even 
those  echoes  composing  the  range  spread.  Al]  the  vertical 
echoes  have  a  positive  doppler  shift  (dark  numbers)  indicating. 


a  downward  motion  of  the  ionization.  The  echoes  between  7  and 
10  MHz  from  the  northeast  have  predominantly  negative  doppler 
frequencies,  the  ionization  is  moving  away.  A  sequence  of 
directional  and  doppler  ionograms  is  required  to  resolve  the 
dynamical  processes  in  the  ionosphere. 
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6.0  REMOTE  TRANSFER  OF  DIGITAL  IONOGRAMS  IN  REAL  TIME 


Remote  display  of  ionograms  has  become  feasible  with 
the  advent  of  digital  ionosondes.  Transmission  of  digital 
ionograms  via  commercial  telephone  lines  poses  no  major  dif¬ 
ficulties  and  allows  to  simultaneously  monitor  ionospheric 
conditions  at  the  various  locations  of  digital  sounders.  The 
need  for  real  time  remote  ionogram  information  arises  in  a 
number  of  circumstances :  monitoring  the  midpoint  of  high- 
frequency  communication  links  or  bistatic  hf  experiments,  mon¬ 
itoring  the  ionospheric  relfection  points  for  over-the-horizon 
radar  operation,  or  simply  the  testing  of  local  field  sites. 
The  distances  over  which  the  transmission  takes  place  may  be 
less  than  100  km  in  the  latter  case  to  more  than  1000  km  in 
the  first  examples.  Simplex  circuit  telephone  lines,  or  al¬ 
ternately  microwave  links,  are  sufficient  for  the  ionogram 
transmission,  even  though  half  duplex  circuits  (or  the  use  of 
reverse  channels)  are  preferable  since  they  allow  handshaking 
procedures  and  data  retransmission  after  error  detection. 

Figure  9  shows  an  ionogram  from  AFGL's  Goose  Bay 
Ionospheric  Observatory,  Canada,  printed  in  Maine,  U.S.A., 
approximately  1000  km  south  of  Goose  Bay  (Smith  et  al,  1979). 
In  the  Centronics  printer  that  was  used  for  the  generation  of 
this  ionogram  the  font-controlling  ROM  (read-only  memory)  was 
replaced  by  a  differently  programmed  chip  in  order  to  produce 
the  optical  weighting  of  the  number  font.  The  time  required 
for  the  remote  transmission  of  an  ionogram  is  g-.ven  by  the 
maximum  baud  rate  of  the  modems  and  the  telephone  lines,  and 
by  the  speed  of  the  printer.  Inexpensive  modems  and  telephone 
lines  limit  the  baud  rate  to  1200  bits  per  second,  or  120 
characters/sec.  For  an  amplitude  ionogram  with  128  range  bins 
(Bibl  and  Reinisch,  1978a)  a  little  more  than  1  sec  is  re¬ 
quired  for  transmission  of  one  data  line.  If  the  printer  has 
a  line  buffer  that  accumulates  the  new  line  while  printing 
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the  previous  one  it  takes  less  than  three  minutes  for  the 
remote  printout  of  an  ionogram  with  160  frequencies  (0  to  16 
MHz  in  100  kHz  increments). 

In  general  it  is  undesirable  to  impose  the  speed 
limitations  of  the  remote  printout  facility  onto  the  sounder 
operation.  We  therefore  store  a  complete  ionogram  in  memory 
and  start  remote  transmission  at  the  end  of  the  ionogram  scan. 
Some  limited  postprocessing  was  performed  in  the  example  shown 
in  Figure  9:  (1)  signals  with  extraordinary  polarization  are 

suppressed,  (2)  all  vertical  echoes  are  underlined,  (3)  digi¬ 
tal  noise  filtering  is  applied.  Of  course  the  specific  type 
of  postprocessing  will  depend  on  the  user's  requirements. 


SUMMARY 


Postprocessing  of  digital  ionograms  has  become  an 
important  step  of  ionospheric  research  with  ionosondes.  The 
new  digital  ionosondes  with  their  multiparameter  data  offer 
new  opportunities  for  research  and  we  have  shown  in  a  number 
of  examples  how  automatic  processing  of  the  data  can  reveal 
specific  ionospheric  features. 
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